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Abstract 



We analyze spin-dependent transport through a spin-diode in the pres- 
ence of spin-flip and under influence of temperature bias. The current po- 
larization and the spin accumulation are investigated in detail by means 
of reduced density matrix. Results show that the spin accumulation is 
linearly increased when the metallic electrode is warmer whereas, its be- 
havior is more complicated when the ferromagnetic lead is warmer. Fur- 
thermore, spin-flip causes that the current polarization becomes not only 
a function of spin-flip rate but also a function of temperature. The current 
polarization is reduced up to 90% if the time of spin-flip is equal to the 
tunneling time. The behavior of spin-dependent current is also studied as 
a function of temperature, spin-flip rate, and polarization. 



1 Introduction 

Study of electron transport through devices fabricated from quantum dots (QDs) 
has attracted a lot of attention during two recent decades [TJ [21 [31 HJ [SJ [HI El 
8,9. Transport through QDs exhibits novel and interesting phenomena such 
as Coulomb and spin blockade effects QUI [TTJ QJ [T3] , Kondo effect [H[T3[T1], 
negative differential conductance [T7J [T5] , and so on. Coupling of QD to differ- 
ent leads such as, normal metal, ferromagnet, or superconductor, is now feasible 
due to recent progress in nanotechnology. One of the most interesting configu- 
rations is a quantum dot coupled to a normal metal(NM) and a ferromagnetic 
(FM) lead. In recent years the configuration has been extensively studied ex- 
perimentally and theoretically 19, 2Q1 [2H [221 123] - With regard to spin-current 
rectification effect observed in the device, it can work as a spin-diode. 

the most studies done about the spin-diode have been focused on the be- 
havior of the system in the presence of electric bias. Due to recent advance in 
the field of thermoelectricity, it is now possible to produce the spin current by 
applying temperature bias across a ferromagnetic semiconductor [241 125] or a 
magnetic insulator [25]. Very recently, F. Qi and co-workers [27] have studied 
the transport through a QD coupled to a normal metal and a ferromagnetic 
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electrode under influence of the temperature bias, they have reported that the 
temperature gradient results in a rectification effect in the current polarization. 
In this article, we analyze the behavior of the system in the presence of spin-flip. 
Transverse magnetic field j^S], spin-orbit interaction and etc can cause the 
spin-flip. With regard to the structure of the device, it is completely possible 
that the electron of the QD interacts with the polarized nucleus of the sublayer. 
Current polarization, spin accumulation, and spin current are analyzed in detail 
by means of reduced density matrix approach [501 IE] • It is found that spin- flip 
affects current polarization significantly. 

In the next section, model Hamiltonian and equations used for calculations 
are presented. Using density matrix approach and wide band limit, formal 
expressions of electron density and current are given. Numerical results and 
analytical expressions for the spin accumulation and the current polarization 
are presented in section 3. In the end, some sentences are given as conclusion. 



AT>0 



AT<0 







Figure 1: Schematic diagram of a spin diode. AT denotes the temperature 
difference between leads. The left lead is warmer if AT > and AT < means 
the right lead is warmer. The increase of temperature causes that the electrons 
of the lead occupy the states above the chemical potential shown as gray space. 



2 Model 

We consider a single level quantum dot coupled to a normal metal (NM) and 
a ferromagnetic (FM) lead as fig. 1. The Hamiltonian describing the system is 
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given as follows 

H = ^ £akacl kcr C akc r + ^ £ <? n <J + U (1) 
aha o 

R[d\di, + d\d t \ + Yy«<™ c ik« d ° + H - c ] 



aka 



where non-interacting quasi-particle approximation is used to describe the elec- 
trodes, c a ka ( c akcr) destroys (creates) an electron with wave vector k, spin a, and 
energy e ak a in the lead a (a = L,R). d a (d <! a ) is annihilation (creation) operator 
in the QD which destroys (creates) an electron with spin a in the QD. U denotes 
on-site Coulomb repulsion and n a — S a d a is occupation operator. The fourth 
term describes spin- flip process and R is spin- flip rate which is spin-independent. 
The last term describes tunneling between the QD and the electrodes and V a ka 
stands for the coupling strength. 

It is clear that the QD can be in : empty state |0 >, singly occupied state 
\a > (e CT ), or doubly occupied state |2 > (S2 = £f + £4. + U). States | f> and 
I 4-> are not the eigenstates of the isolated QD Hamiltonian because of spin-flip. 
Using Markov approximation the time evolution of the density matrix elements 
are given as [301 OH] 

= -i < s\[H s . f ,P]\s' > +6 SS , Y} w ksPkk - W sk P ss ] (2) 

k^s 



dP ss > 
dt 



hi- s ss ,\[J2 w sk + W s'k\P S S> 



k^s s'^k 

where H s j denotes the fourth term in Eq.((T]). P ss is probability of being in 
the state s (s = 0, a, 2) whereas, P ss > (s ^ s') describes coherency between 
states I t> and | J.>. W ss > stands for transition from state \s > to \s' > and is 
computed by Fermi's golden rule as 

W aaf =J2^[fa(\e a ,'\)S N .,N^+l + f a {\E^\)S r f.,Ni-l] (3) 



fa(x) — (1 + exp((x — Ma)/^s^a)) _1 is Fermi distribution function of lead a 
where T a and \x a are temperature and chemical potential of lead a and f a = 
1 — f a . is spin-dependent tunneling rate of the lead a obtained by using 
wide band approximation. As an instance, Pfi is obtained from Eq.([2|) as 

^ti = iR [P n P U ] - i[£ W n + W lk ]P n (4) 

and Ptf= p n- 

Now, we define n a = P aa + P 2 2- Using Eq.Q and normalization condition 
(-P00 + So- P<t<t + P22 = 1), it is straightforward to show 

— f = W a [1 - n a - n g ] - W aQ na + W s2 n s + iR[P aS - P Sa ] (5) 
dt 
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where a is opposite of a. Solving Eq.([S]) in the steady state = 0), the 

spin-dependent current crossing from the lead a is given as 

% = W£ a [l-n a - n B ] - W^n a + W« 2 n s (6) 

In the continue, we consider r = IQpeV as lead-QD tunneling rate and 
set = T and T% = F [l + ap] where a = 1(-1) for t (|), and p is the 
spin polarization degree of the ferromagnetic lead. Note that the left lead is a 
normal meta thus its tunneling rate is spin- independent whereas, the right lead 
is a ferromagnet whose majority carriers are assumed to be spin-up electrons. 
In addition, we assume e a = -2meV, U = 5meV, and AT = T L - T R . AT > 
means the left lead is warmer (Tl = Tq + AT) while, AT < means Tr = 
To + AT where Tq = 0.2meV is the base temperature of the system. 




Figure 2: Spin accumulation versus temperature bias. 5 = and p = 0.4 (solid), 
S = and p = 0.8 (dashed), S = 0.5 and p = 0.4 (dotted) and S = 0.5 and 
p = 0.8 (dash-dotted). Inset shows a (solid), /3 (dashed), and a + (3 (dotted). 



3 Results and discussions 

In order to analyze the spin accumulation and the current, we have assumed 
that f R {e a ) = 1, f R {e a + U) = 0, f L (e a ) = a and f L (e a + U) = p for AT > 
whereas, fh{sa) = 1, !h{e a + U) = 0, /i?(e CT ) = a and f R (e a + U) = ft in 
AT < 0. Note that a and (3 are temperature-dependent and their dependence on 
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the temperature is shown in inset of fig. 2. Fig. 2 shows the spin accumulation 
(m = n-f — nj.) as a function of temperature bias. Results show that the most 
change in m is happened when the temperature difference is very low. As one 
can see in inset of fig. 2, a — > 1 and /3 — > in low temperature, so that there is 
only one electron in the QD and as a result, the ferromagnetic electrode affects 
the spin properties of the system significantly. By increase of temperature, the 
second electron enters the QD thus the spin characteristics of the system are 
vanished. By solving Eq.© in the steady state we obtain blow equation for the 
spin accumulation as a function of temperature, polarization, and spin-flip rate. 



2p[(a - l) 2 - /3 2 



4S 2 [3 + a-P]+y[4-(a+l- /3) 2 ] 
2p[a + P-l]y 



AT > (7a) 
AT < (7b) 



P 2 y[(a - p) 2 - 1] + AS 2 [Z + a -0\- y[(a - /3 + l) 2 - 4] 

where y = 1 + j3 — a. 

From Eq.©, it is obvious that increase of temperature difference results 
in m = because of a, /3 —¥ 1/2. Physically, by warming an electrode, elec- 
trons near the Fermi surface occupy states above the chemical potential, see 
fig. 1, so that there will be some holes below the chemical potential. The 
excited electrons have the necessary energy to overcome the charging energy 
so that the QD will have two electrons and as a consequence, m = 0. Eq.flT]) 
also shows that the dependence of m on p and S is different in positive and 
negative temperature bias and this difference results in an asymmetry in spin 
accumulation-temperature bias characteristic of the system. As one can see in 
Eq.(j7|), m decreases rapidly by increase of S hence the variations of the spin 
accumulation become inconsequential in the presence of spin-flip, as it is shown 
in fig. 2. In positive and low temperature difference, m is positive because the 
ferromagnetic lead acts as an emitter and due to T^ R > T^, the probability of 
being in the spin- up state is more. Note that a + (3 < 1, see dotted line in the 
inset, and as a result (a — l) 2 — (3 2 becomes positive in Eq.flT]). In negative and 
low temperature difference, the left lead acts as the emitter and as a result, m 
becomes negative because the spin-down electron has to stay in the QD for a 
longer time. 

The dependence of the spin accumulation on the spin-flip rate and the po- 
larization is plotted in fig. 3. asymmetry of m versus R is well seen in fig. 3a. 
As expect from Eq.Q, rn is significantly reduced by increase of R. The depen- 
dence of m on p is much more interesting. In positive temperature bias, the 
spin accumulation is linearly increased by increase of p whereas, its behavior is 
completely different in negative temperature difference. The main result is risen 
due to the change of the emitter. In AT > where the right lead is emitter, the 
linear relation is dominant while, in the condition that the left lead is emitter, 
the dependence is more complicated. We reported the same behavior for the 
spin accumulation in a spin-diode under influence of voltage bias |32j . 

Spin-resolved currents as a function of temperature bias, spin-flip rate, and 
polarization are shown in figs. 4 and 5, respectively. First, we estimate the 
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current in the presence of positive temperature bias. By solving Eq.Q and 
with respect to Eq.(|6]), spin-dependent currents are obtained as: 

l\ = T [a-r H -(a-p)n i ] (8a) 
li=T [a-(a-j3)n t - ni ] (8b) 

where n a is given as 

pa[(l - a) 2 - /3 2 ] + 4S 2 (1 + a) + (g + l) 2 + q[(a - /3) 2 - (1 + a)] 
n<T 4S 2 [3 + a-/3] +y[4- (l + a-/3) 2 ] lj 

where a = 1(— 1) for j" (4.). In positive bias, the current crossing from the left 
lead is negative because the right lead acts as the emitter, as explained before. 
From Eq.©, it is clear that n-j- is reduced by increase of R whereas, increase of 
R gives rise to increasing n^. This effect results in reduction of 1^ and increase 
of Jj_, as it is seen in Eq.©. Note that the probability of finding the QD in 
spin- up state is more in positive temperature bias if spin- flip process is absent. 
In the presence of spin-flip, the spin-up electron may rotate and change to a 
spin-down electron. Therefore, the spin-flip increases the probability of being in 
the spin-down while, the probability of being in spin-up state is reduced. With 
respect to Eq.©, although the existence of R influences significantly on the 
current, the magnitude of it does not have significant role in the curvature of 
the current due to existence of a large term in denominator. Furthermore, 1^ 
will be equal to 1^ if R » Tq because of rif = n^, see Eq.©. From Eqs.(8,9), 
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it is obvious that the current is related to p linearly. This dependence is well 
observed in fig. 5b. It is very interesting to note that unlike R, increase of 
p results in increase of V and decrease of -Tk Indeed, with increase of p the 
spin- up electron is injected into the QD faster and as a result U > 1^. 

Now, we analyze the current when the ferromagnetic lead is warmer. In this 
case, is equal to 1^ and given as 

r =r„[i-»t-»4.] ( 10 ) 

where n a is obtained from 

p<ry[a + (3 - 1] - ay 2 p 2 + 45 2 [1 + a] + y 2 [l + a] 
Ua p 2 y[(a-p) 2 -l}+4S 2 ('3 + a-P)-y[(l + a-f3) 2 -4} 1 ' 

The first result observed in fig. 4 is the current is positive i.e. the left lead 
acts as the emitter. It comes from the fact that electrons below the Fermi level 
in the right lead are lesser than electrons in the left lead. On the other hand, 
increase of R results in increase of both F and 1^. From inset of fig. 4, it is 
clear that increase of R gives rise to reduction of and increase of n-j-. Note 
that the probability of being in the spin-down state is more if R — because 
the spin-up electron injected from the left lead leaves the QD faster than the 
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Figure 5: (a) I a as a function of spin-flip rate. We set p = 0.8 and AT = ±3. 
Solid line is V and dashed line is J-k AT = —3 is plotted in gray, (b) spin- 
dependent current versus p. We set R = 0. Inset shows rif (solid) and nj, 
(dashed) for AT = -3. 



spin-down electron. As it is obvious in inset and Eq. (|ll[) , the rate of reducing 
n± is faster than the rate of increase of rif. This fact leads to increase of I" in 
the presence of spin-flip. Increase of I a due to spin-flip is clearly seen in fig. 
5a. Like AT > 0, the spin-resolved currents saturate for R > 0.5Fo. Unlike 
AT > 0, the dependence of I" on the polarization is not linear anymore. This 
dependence is plotted in fig. 5b. The behavior of I a is very interesting in p = 1. 
In positive temperature bias, J* becomes zero because there are no spin-down 
electrons in ferromagnetic lead to enter the QD, but V = I* = if AT < 0. It 
is straightforward to show that nf + nj. = 1 if S = and p = 1 and with respect 
to Eq. tjTU)) . I a becomes zero. Indeed, the electron inside the QD cannot tunnel 
to the right lead so, other electron cannot enter the QD. In addition, increase 
of p results in decrease of I a . 

The current polarization, £ = ffifjz M plotted in fig. 6a. With respect to 
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Figure 6: Current polarization versus (a) temperature bias, (b) polarization, 
and (c) spin-flip rate. 



Eqs.(8,10), it is straightforward to show that 

C=^H^ f° r AT > Q ( 12a ) 
C = for AT<0 (12b) 

Note that our results are identical results given in Ref. [27] if S = and a = 
/3 = 1/2. Presence of spin-flip causes that the current polarization becomes not 
only a function of spin-flip rate but also a function of temperature because of 
presence of a and (3 in Eq. (fT2"]) . As expected, spin- flip gives rise to reduction of 
£ because this process tries to destroy the spin characteristics of the system. If 
spin-flip rate becomes equal to tunneling rate, ( is reduced up to 90%. However, 
the system can still work as a rectifier although its performance is weak. The 
dependence of C on polarization and R is shown in figs. 6b and 6c, respectively. 
Linear behavior of £ versus p is well seen. £ = 1 if R — and p = 1 because 
there is no spin-down current through the system but, £ is significantly reduced 
in the presence of spin-flip because the spin-up electron into the QD can change 
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to a spin-down electron and as a result, the spin-down current is created. Unlike 
fig. 6b, the behavior of C against R is nonlinear. Such behavior about current 
polarization was before reported in the presence of electric bias [32] . 

4 Conclusion 

In this article, We study spin-dependent transport through a quantum dot cou- 
pled to a normal metal and a ferromagnetic electrode. Analytical expressions 
for spin accumulation and current polarization are obtained in the presence 
of spin-flip, using master equations. Results show that the dependence of the 
system on the polarization is linear when the metallic lead is warmer whereas, 
its behavior is more complicated when ferromagnet is warmer. On the other 
hand, the spin-flip affects significantly the current polarization so that it is re- 
duced up to 90% if the spin-flip time is equal to tunneling time. The effects of 
temperature, polarization, and spin-flip rate on the spin-resolved currents are 
also estimated. When the metallic electrode is warmer, spin-flip results in the 
increase of spin-down current while, spin-up current is decreased. 
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